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PART  I 


IR  RADIATION  MODEL  FOR  EXCITED  OZONE 
PRODUCED  IN  A PURE  02  ELECTRIC  DISCHARGE 

Kurt  L.  Wray  and  George  E.  Caledonia 


Physical  Sciences,  Inc. 

Wakefield,  Massachusetts 

Introduction 

The  theoretical  work  reported  here  was  done  in  conjunction  with  the  experi- 
mental program  being  carried  out  at  the  Astronomy  Research  Facility  at  the 
University  of  Massachusetts  under  the  direction  of  Prof.  John  Strong.  The 
experimental  program,  directed  at  the  study  of  the  chemistry  of  species  produced 
by  electrical  discharges  in  a simulated  upper  atmosphere,  is  being  carried  out 
in  a 100  feet  long,  3 feet  diameter  discharge  tube.  The  tube  has  cylindrical 
coaxial  geometry  with  the  outer  wall  serving  as  one  electrode  and  a 33  feet  long, 

3 inch  diameter  tube  located  axially  along  the  central  one  third  of  the  larger  tube 
serving  as  the  other  electrode.  The  optical  system  allows  triple  pass  operation 
and  the  test  gas  can  be  studied  in  either  absorption  or  emission.  Neglecting 
axial  diffusion,  the  effective  optical  path  for  discharge  produced  species  is  then 
100  feet  or  3 x 10^  cm. 

Experimental  Conditions 

A visibly  uniform  discharge  was  achieved  with  600  volts  (direct  current ) at 
a pressure  of  0.1  Torr;  the  total  current,  I,  drawn  under  these  conditions  was 
0.3  amp. 

The  electric  field,  E,  for  cylindrical  geometry  is  given  by 


E » 


r In  (rg/i’j)  2 . 5r 


(1) 


where  V = 600  is  the  potential  across  the  electrodes  and  ^ = 3.8  cm  and 

r2  = 46  cm,  the  inner  and  outer  electrode  radii,  respectively.  At  a pressure  of 

0.1  Torr, 


57  < E/P  < 720  (volt  cm"1  Torr"1)  . 


(2) 


In  pure  02  this  corresponds  to  an  electron  temperature  Tg  of 


3 


T > 2 . 5 eV  . 
e ~ 

If  there  were  no  ionization  in  the  gas,  the  electron  concentration  ne  would  be 


(3) 


given  by  the  equation 


ne  " 2 7ir  1 v Q 


(4) 


where  1 = 1 x 103  cm  is  the  discharge  length  and  Q is  the  charge  on  an  electron^ 
6 x 10“19  coulomb) . The  electron  drift  velocity  v was  taken  to  be  v - 2x10 
cm  /sec  as  typical  for  the  discharge  conditions,  hence,  for  I = 0.3  amp 


S 3 

n « 3x  10  electrons/ cm 
e 


(5) 


at  the  outer  wall.  It  is  probably  not  much  more  than  a factor  of  10  larger  than 
that  anywhere  in  the  discharge  volume  even  though  the  inclusion  of  gas  ionization 
effects  would  modify  the  r_1  dependence  of  Equation  (4) . In  making  numerical 
calculations  below,  we  arbitrarily  use 


7 3 

n -2x  10  electrons  cm 
e 


< 6) 


as  an  upper  bound  to  the  local  electron  concentration  in  the  present  experimental 
apparatus . 

Heterogeneous  reactions  occurring  on  the  walls  of  the  apparatus  are 
important.  These  are  kinetically  limited  by  diffusion  of  the  reactive  (neutral) 
sp  cies  to  the  walls.  For  the  system  geometry,  the  characteristic  diffusion 
time  rdiff  is  given  by  a complex  transcendental  equation  involving  Bessel 
functions:  for  the  actual  ^ and  r2  involved,  we  find 


1.5  P (sec) 


(7) 


where  P is  the  pressure  in  Torr. 

Kinetic  Mechanisms 

Table  I gives  the  kinetic  processes  to  be  considered  in  a pure  C>2  discharge 
for  the  experimental  conditions  discussed  above.  Also  given  in  Table  1 are  the 
rate  constants  k for  2-  and  3-body  reactions  1 through  6,  or  time  constants  r for 
unimolecular  processes  7 through  9 . Here  we  comment  on  each  of  the  rate 
constants  given  in  Table  1 . 
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Reaction  1:  Discussed  below . 

2:  Well  known. 

3:  Recently  measured  (C.  von  Rosenberg). 

4:  Unknown  but  process  probably  nut  important  here  because  both 

reactants  are  trace  species. 

5:  Well  known. 

6:  Well  known. 

7:  Diffusion  calculation  is  accurate  and  y » 1 is  good  assumption, 

8:  y is  unknown  but  process  is  unimportant  at  P > .03  Torr  even 

with  ysl. 

, . 'system 

9:  Assumes  tr  (v  = 4 - 3) ^ 

Reaction  1 is  listed  in  Table  I as  a global  reaction.  The  basis  for  this  can 
ie  seen  in  Figure  1 where  the  rate  constants  for  the  following  reactions  are 
dotted  vs.  E / P. 

Reaction  10:  e + C>2  — O + O 

11:  O"  + 02  - O + 02  + e 

12 : e + 02  - 02  + + 2e 


1 : e + Or 


TABLE  1 . PURE  02  DISCHARGE  KINETICS 
Global 


2 : O + 2 02  — 
3:  03  + 02  - 

4:  03  + O - 
5:  03  + O - 
6:  O+O  + O, 


20  + e 


O3  + o2 

°3  +02 


2 02 
2 02 
2 Or 


3 -1 

k = f (E/P)  cm  sec 

k = 6 x 10"34  cm6  sec-1 

-14  3 -1 

k = 2 x 10  cm  sec 


k = 6 x 10"15  cm3  sec-1 
k = 5 x 10"33  cm6  sec"1 


a =1 

7:  O 1/2  09 

Wall  2 

8:  °3  Waif  °3 


9:  Or 


O3  + hy 


td  = 1.5  PTorr  sec  (system) 


rD  = ll5a' 


PTorrSeC  (a<  1} 

v-2 


R 


= 2.5  x 10"Z  sec  (v  = 4 - 3) 


* Vibration  ally  excited. 
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It  ic  noted  that  Reaction  10  is  rate  limiting;  for  the  E/ P range  of  interest,  the 
0“  produced  is  rapidly  destroyed  by  collisional  detachment,  the  sum  of  11  and  12 
being  the  global  reaction  1 with  an  effective  rate  constant  of  ~2  x 10"11 . 

Reaction  12  identifies  02  as  the  dominant  positive  ion,  whose  concentration  is 
equal  to  the  electron  concentration  since  negative  ions  are  unimportant  here. 

We  have  not  included  in  our  model  ozone  destruction  by  electrons: 


Reaction  13:  e + O3 


■•O  + 0«l 


‘ o~2  + o 


-19  2 

which  processes  have  cross  sections  cj  ~ 10  cm  at  electron  energies  ~ 5eV. 
This  cross  section  would  make  Reaction  13  unimportant  for  ne  £ 10  cm 

Furthermore,  we  have  not  included  in  our  model  direct  (direct  or  indirect) 
electron  excitation  of  vibration: 

Reaction  14:  e + 03  - e + 03  k = unknown 

15a:  e + 02  - e + 02 

15b:  02  +03-63  + 02 


These  reactions  are  not  considered  to  be  important- -more  will  be  said  about 
Reaction  14  below. 

Steady  State  Solutions 

The  rate  equations  for  [ 0|  , [ O3  ] . and  | O3]  are  given  in  Table  2.  Here 
[ ] are  mole  fractions  and  N is  the  total  number  of  particles/ cm^  which  is 
essentially  the  02  concentration.  At  the  bottom  of  Table  2,  complete  equations 
are  given  for  the  steady  state  mole  fractions  obtained  by  setting  d[  ] / dt  s 0. 
Crossed  out  terms  were  shown  to  be  negligible  for  present  conditions.  The 

if 

radiative  and  catalytic  decay  of  O3  are  readily  shown  to  be  unimportant  at  the 
higher  pressures;  i.e. , P it  0.06  Torr  for  radiative  decay  and  P 0.03  Torr 
for  wall  catalysis.  Hence  at  P > 0.06  Torr,  assuming 


[Ol  « 1 o2| 

and  using  the  rate  constants  of  Table  1,  we  find  (to  1 significant  figure) 

3 P k.  nn  , . ..  „ 

f ol  = 1 g catalytic  wall 

1 1 ' . _3  for  O-atoms 

1 + a P 


(8) 

(9) 


* Vibrationally  excited. 
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TABLE  2 . KINETIC  EQUATIONS  AND  STEADY 
STATE  SOLUTIONS 

d[Oj  = 2^  - k2[  01  N2  - k4  [O3]  [ O]  N - k5[03]  fOl  N 

- 2kg [ 01  2N2  - !£] 
tD 


= k2[  01  N^  - k3 


[03*]N-k4[0;]|01N-^.I^] 


r *-1  r -1  o\o*]  [0*  ] 

= k3  [°3  j N - k5  [°3]  I °]  N + 4^  + 4^ 


Whence  in  steady  state 


[ O) 

K] 

hi 


^N2  +k4[o']N  +k5[03]N+2kei^tl3r+TD 


k2[  0|  N2 


k5|OlN 


* Vibrationally  excited. 
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where  P is  the  pressure  in  Torr  and  a varies  monotonically  with  time  ircm 
1 to  2, 


[03I 

= 1 x 10"3  P [ O]  , 

(10) 

[03] 

= 3.3  x 10"3  P . 

(ID 

In  Equation  (9)  the  1 in  the  denominator  results  from  the  recombination  of 
O-atoms  on  the  walls  of  the  discharge  tube.  For  non-catalytic  walls  (y  ~ 0 for 
Reaction  7),  Equation  (9)  becomes 

3 k . n 

[O]  = — (non-catalytic  wall)  . (12) 

a P2 

1 1 7 

Thus,  for  P = 0.1  Torr,  k.  = 2 x 10  and  arbitrarily  taking  n = 2 x 10 
o 1 e 

cm  we  find  for  catalytic  walls 

[ O]  = 1.2  x 10“4 

[O*  | = 1.2  x 10"8  (13) 

[03]  = 3.3  x 10-4  . 

Note  that  [ O3  ] / [ O3]  = 4 x 10"5,  which  is  considerably  smaller  than  the 
thermal  equilibrium  population  at  300°  K,  where  [ O3  ] /[  O3]  = 6.7  x 10  3. 

For  non-catalytic  walls, 

[ O]  = 1 .2  x 10"1 

[O3  ] = 1.2  x 10"5  (14) 

[ O3]  = 3.3  x 10"4  . 

The  steady  state  mole  fraction  of  ozone  [ O3]  is  independent  of  [ O]  (and 
hence  wall  catalycity)  because  O3  is  made  by  O-atoms  (Reactions  2 and  3) , but  it 
is  also  destroyed  by  O-atoms  (Reaction  5) . On  the  other  hand,  O3  is  made  by 
O-atcms  (Reaction  2)  but  it  is  destroyed  by  Og-molecules  (Reaction  3)  and  hence 
[ O3]  scales  like  [ O]  . 

* Vibrationally  excited. 
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Time  History 


Approximate  time  histories  for  [ O]  , [Og]  and  [ O3  ] are  shown  for  the 
catalytic  case  in  Figure  2 . Here  the  slopes  (d  [ X]  /dt)  are  established  by  the 
dominant  term  (s)  in  the  differential  equations  of  Table  2,  whereas  steady  state 
levels  are  approached  at  about  3 destruction  times  given  by  the  dominant  negative 
term  (s) . It  should  be  noted  that  O reaches  steady  state  in  about  0.4  sec  but  Og 
takes  15  minutes. 

Radiation  Predictions 

■ _3 

The  IR  radiation  emitted  in  the  instrument  band  pass  ( ^ Ag  p =10  /jl  ) 
by  the  Og  can  be  approximately  calculated  from  the  equation 


he 

A 


AAB.P. 
AA  Band 


(15) 


For  the  experimental  conditions  with  catalytic  walls,  we  have: 


. _ 1.2  x 10'8  x 3.5  x 1015  , he  .3xl03xl£^  (16) 

03  4 tt  x 0 . 1 A 1 

-12  2 

— 2x10  watts  cm  - ster. 

3 

For  non-catalytic  walls,  Iq*  would  be  10  times  larger  . 

These  intensities  must  be  compared  with  the  threshold  sensitivity  of  the  IR 
optical  system  which  is  part  of  the  experimental  facility 


*th 


reshold 


4 x 


10"8  watts 


cm  ster 


(17) 


Hence,  it  is  seen  that  the  present  system  misses  seeing  the  discharge  produced 
IR  signal  by  about  4 orders  of  magnitude.  Even  if  the  walls  of  the  discharge 
tube  were  made  non-catalytic  to  O-atoms  the  system  would  still  be  inadequate. 
We  believe  our  calculations  of  [ 03  | to  be  an  upper  bound. 

It  should  be  pointed  out  tlat  nothing  much  can  be  gained  by  varying  the 
pressure.  Note  that  for  the  non-catalytic  wall  case  (Equations  (10)  and  (12)) 


3 x 10"3  kj  ng 
a P 


for  [ O]  << 


(18) 


* Vibrationally  excited. 
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(19) 


so  that 

■o’  “ [03]  N*ne  • 

Hence  it  seems  that  the  signal  could  only  be  increased  by  increasing  ne,  but  note 
that  Equation  (12)  gives 

[O]  ocHl  . (20) 


which,  for  an  assumed  ne  of  2 x 107  was  already  about  as  large  as  allowable 
( [ O]  > 0.1 ) . We  therefore  conclude  that  a signal  gain  vould  only  be  accom- 
plished by  simultaneously  increasing  ng  and  P such  that  [ O]  at  constant. 

This  is  probably  not  possible  to  do. 

* 

e + Og  — e + Og 

As  already  mentioned  above,  the  rate  constant  for  Reaction  14  is  not  known. 

For  the  conditions  considered  above,  however,  it  is  readily  shown  that  Reaction 

* 

14  would  dominate  the  Og  production  only  if 

k14  ~ 10  kl  (catalytic  wall)  (21) 

k14  > 104  k1  (non-catalytic  wall)  . (22) 

It  is  most  unlikely  that  k14  Z.  10"9,  hence  Equation  (22)  can  not  be  achieved  and 
Reaction  14  is  unlikely  to  affect  the  conclusions  reached  above. 

Finally,  let  us  ask  if  the  present  experiment  could  be  modified  to  measure 
the  rate  of  Reaction  14.  Suppose  pure  Og  were  used  as  the  test  gas  in  the  dis- 
charge tube,  then,  an  upper  bound  to  steady  state  [ Og]  is  given  by 


k14ne 
k3  N 


(23) 


where  K = k,  = 2 x 10  14,  P = 0.1  Torr  of  0«,  nc  s 2 x 107  (upper  bound) , and 
3 o q o e 

k14  < 1 x 10  y,  we  find  that 

[o*]  < 3x  10-4  . (24) 


* Vibrationally  excited. 
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This  [ o!  1 is  about  4 orders  of  magnitude  larger  than  that  predicted  for  the  case 
where  ol  is  the  test  gas.  However,  that  case  produced  insufficient  O3  by  about  a 
factor  of  104.  Considering  that  Equation  (24)  represents  an  unlikely  upper 
bound,  the  0~  test  gas  experiment  is  not  considered  likely  to  succeed.  In  any 
case,  the  [ O*  1 concentration  shown  in  Equation  (24)  would  be  less  than  room 

temperature  thermal  concentration. 


* Vibrationally  excited. 


13 


PART  II 


EXPERIMENTAL  PROGRAM 

J.  Strong  and  P.  Hansen 

Astronomy  Research  Facility 
University  of  Massachusetts 
Amherst,  Massachusetts 

Experimental  Facility 

The  discharge  column  was  contained  in  a large  steel  absorption  cell  1 meter 
in  diameter  and  30  meters  long.  The  inside  of  the  cell  is  normally  coated  with 
a nonconductive  epoxy  paint.  Thus,  a 12 -meter  section  in  the  center  of  the  cell 
was  coated  with  conductive  paint.  This  formed  the  outer  electrode.  The  inner 
electrode  was  a 7.5  cm  diameter  12  meters  long  aluminum  tube  centrally 
suspended  in  the  cell.  By  triple  passing  the  cell,  the  effective  optical  path- 
length  (Figure  3)  of  the  discharge  column  was  36  meters. 

The  discharge  was  powered  by  a 5 Kw  electronically  controlled  power 
supply.  Modes  of  operation  included  either  d.c.  or  60  Hz  a. c.  Good  steady 
discharges  were  obtained  in  the  0.05  to  0.35  Torr  pressure  range.  Alternating 
current  discharges  were  more  stable  than  d.c.  discharges  due  to  arc  suppres- 
sion. A Langmuir  probe  was  used  to  sample  electron  temperatures  when  the 
discharge  was  operated  in  the  d.c.  mode.  Electron  temperatures  in  all  cases 
where  a discharge  could  be  maintained  were  found  to  lie  in  the  range  of  0.5  to 
2.5  ev.  A residual  gas  analyzer  was  used  to  monitor  impurities  in  the  discharge 

The  spectroscopic  equipment  included  a Hadamard  transform  spectrometer 
(Reference  1)  (used  in  the  single  scan  mode)  covering  the  wavelength  range 
from  6 to  21  microns,  a Brower  M10  monochromator,  and  a Perkin-Elmer 
210B  monochromator.  The  latter  two  instruments  were  used  to  cover  the  1 to  4 
micron  wavelength  range.  Both  monochromators  were  modified  for  use  with 
cooled  PbS  detectors. 

Measurements 

* 

Our  emphasis  has  been  on  attempts  to  observe  excited  ozone  Og  in  air  and 
oxygen  discharges.  The  discharge  was  excited  with  both  direct  and  alternating 


1.  P.  Hansen  and  J.  Strong,  Appl.  Opt.  11^,  502  (March  1972). 
* Vibrationally  excited. 


14 


th  in  the  30-meter  absorption  cell 


current  between  the  outer  wall  and  a central  electrode.  This  discharge  con- 
figuration required  a conductive  outer  wall  which  was  catalytic  such  that,  as 
pointed  out  in  the  analysis  by  Wray  (appended),  it  reduced  the  03  concentration 
by  an  estimated  three  orders  of  magnitude  (as  compared  to  glass  walls) . 

Background  radiation  causes  a serious  problem  because  the  optics  of  the 
discharge  column  are  at  room  temperature  ~ 290K.  This  is  true  especially  for 
the  9 .6p  ozone  band,  a wavelength  where  the  290K  blackbody  curve  peaks.  The 
O3  concentration  could  be  increased  by  increasing  number  density  of  electrons 
and  gas  pressure,  simultaneously.  With  the  present  catalytic  walls  and  discharge 
configuration,  this  could  be  done  by  employing  a thermionic  central  electrode. 

The  resulting  gain  in  the  03  concentration  is  not,  however,  judged  enough  to 
bring  success--at  most,  one  order  of  magnitude. 

Other  methods  of  exciting  discharges  for  the  ozone  studies  have  been  con- 
sidered. One  would  be  to  make  the  outside  wall  non-catalytie . Others  include 
RF  and  spark  gap  discharges  in  order  to  enhance  very  greatly  number  densities 
at  higher  pressures.  RF  discharges  did  rot  appear  feasible  for  us.  Spark  and 
continuous  discharges  using  low  surface  area  auxiliary  internal  electrodes  appear 
to  have  some  merit,  but  would  have  required  considerable  alterations  of  the 
present  discharge  configuration. 

The  best  use  of  the  present  system  lies  in  the  wavelength  range  below  4p 
v'here  background  radiation  is  less  pronounced.  To  explore  here,  we  have  taken 
spectra  of  oxygen,  nitrogen,  and  air  discharges  using  a circularly  variable 
filter  for  low  resolution  spectra,  and  the  grating  monochromators  for  better 
resolution  ( An  ~ 5 cm-1 ) . Pressures  ranged  from  .05  to  .35  Torr. 

Typical  monochromator  spectra  are  shown  in  Figure  2 . The  spectra  were 
obtained  with  a.c.  discharges  with  total  currents  near  3 amperes  and  pressure 
near  .25  Torr.  These  conditions  optimized  the  signal-to-noise  ratio.  Higher 
pressures  caused  non- radiative  dark  spaces  in  the  discharge  column  reducing 
the  effective  optical  cross  section.  At  lower  pressures  p ^ 0.1  Torr.  necessary 
current  densities  could  not  be  obtained  to  generate  appreciable  signals. 

The  oxygen  spectrum  shows  the  atomic  oxygen  lines  previously  seen  by 
Saum  and  Benesch  ( Reference  2 ) . Line  4 in  the  spectrum  is  a previously 


2.  K.A.  Saum  and  W.M.  Benesch,  Appl.  Opt.  9,  1419  (June  1970) . 
* Vibrationally  excited. 
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unobserved  atomic  oxygen  line.  The  nitrogen  spectra  show  the  presence  of  CO 
with  some  as  yet  unidentified  species  near  2.10  and  2.45  microns.  The  air 
spectra  are  similar  to  the  nitrogen  spectra  with  less  s/ n. 

Near  the  termination  of  this  contract,  a cooperative  effort  was  undertaken 
with  R.  Murphy  of  AFCRL,  Bedford,  Massachusettes,  to  obtain  spectra  of  the 
12 -meter  discharge  column  in  the  1 to  6 micron  wavelength  region  using  an 
interferometer  and  cooled  PbS  and  PbSe  detectors.  Background  signal  problems 
were  reduced  by  using  a pulsed  discharge  on  the  chopper. 


17 


PART  m 


EPILOGUE* 

The  effort  described  in  the  preceding  portions  of  this  report  was  originally 
conceived  and  begun  in  order  to  rapidly  acquire  certain  information  concerning 
the  properties  of  irradiated  air-like  gaseous  mixtures,  using  an  already  existing, 
operable  laboratory  device. 

Unfortunately,  this  effort  was  mitigated  by  time,  circumstances,  and  other 
more  productive  programs.  Consequently,  this  laboratory  endeavor  was  termi- 
nated in  early  calendar  1974  by  mutual  but  reluctant  agreement  of  all  parties 
involved. 

The  efforts  of  all  personnel  on  this  project  were  most  commendable  through- 
out its  brief  ev;stence. 


♦Added  in  press  by  the  Defense  Nuclear  Agency's  Contracting  Officers 
Repiesentative. 
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